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Abstract 

Recollapsing homogeneous and isotropic models present one of the key ingredi- 
ents for cyclic scenarios. This is considered here within a quantum cosmological 
framework in presence of a free scalar field with, in turn, a negative cosmological 
constant and spatial curvature. Effective equations shed light on the quantum dy- 
namics around a recollapsing phase and the evolution of state parameters such as 
fluctuations and correlations through such a turn around. In the models considered 
here, the squeezing of an initial state is found to be strictly monotonic in time during 
the expansion, turn around and contraction phases. The presence of such mono- 
tonicity is of potential importance in relation to a long standing intensive debate 
concerning the (a) symmetry between the expanding and contracting phases in a rec- 
ollapsing universe. Furthermore, together with recent analogous results concerning 
a bounce one can extend this monotonicity throughout an entire cycle. This pro- 
vides a strong motivation for employing the degree of squeezing as an alternative 
measure of (quantum) entropy. It may also serve as a new concept of emergent time 
described by a variable without classical analog. The evolution of the squeezing in 
emergent oscillating scenarios can in principle provide constraints on the viability of 
such models. 



1 Introduction 

Classical cosmology has shown enormous progress over the recent years. Despite this a 
number of fundamental questions remain. Central among these is the fact that within the 
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classical general relativistic framework, the initial state of the universe is singular which 
would result in the breakdown of laws of physics. To obtain a satisfactory scenario with 
a non-singular initial state one often looks to quantum gravity and quantum cosmology. 
In fact, with a loop quantization one can generically resolve the big bang singularity in 
cosmological and other models [HE]. In the simplest cases, a bounce results which keeps 
the volume non-zero and the universe away from the classical singularity reached otherwise 
at the big bang. The possibility of a contracting phase (or several phases) before the 
hot big bang has recently been invoked in a number of cosmological scenarios, including 
several models proposed as alternatives to standard inflation, such as for example pre-big 
bang [3] and the ekpyrotic/cyclic scenarios [UEEE]. The assumed nature of such phases, 
however, has so far been mostly rather ad hoc, without a satisfactory treatment of the 
classical singularity. The presence of such phase (s) raises important questions, including 
their nature and their relation to the present phase of the universe. This in turn relates to 
fundamental questions such as, among others, cosmological entropy and the arrow of time. 

Now given that the big bang was a high-energy, strong-curvature regime, the under- 
standing of the pre- and post-bounce phases would require a full control of dynamical 
evolution of the quantum state through such a bounce. Moving through a bounce, a wave 
packet can spread and deform significantly, implying that the universe before the bounce 
could, for all we know, have been in a state very different from what we see now. Thus, 
to understand the cosmological dynamics through such bounces, all aspects of a quantum 
space-time are essential, including its fluctuations and higher moments. 

In loop quantum cosmology, solvable models with controlled state properties exist if 
the matter source is a free, massless scalar. This has been analyzed numerically [7J [H [9j 
[TO] and analytically [TTI [T2] . More general models can be treated by means of effective 
equations [131 E], as they are also employed here for the recollapse. Note that the concept 
of effective equations is much more general than simply providing correction terms to 
classical equations. With a complete set of consistent effective equations one can, in fact, 
derive dynamical properties such as expectation values, fluctuations, correlations or higher 
moments for full quantum states. As we will see below, state properties can be studied 
directly by using effective equations, which provide an economical and representation- 
independent approximation scheme of the evolution of states. (For another discussion of 
effective equations especially in quantum cosmology, see [15].) 

If one combines the quantum bounce with a classical recollapse, cyclic models ensue. 
Such oscillatory models, according to which the universe undergoes many (and possibly 
an infinite number of) bounces, have been employed in order to construct non-singular 
emergent models which can set the initial conditions for a successful phase of inflation. 
Since such a universe can pass through many cycles, and hence many high energy, strong- 
curvature regimes, this could result in even more severe changes of its state compared to a 
single bounce. We should note that oscillatory models have a long history in cosmology at 
least since the studies by Tolman in the 1930's [TH [T7J - albeit within a classical setting. 
Interestingly, Tolman also considered the question of cosmological entropy for these models, 
claiming that the entropy during the expanding phase should be slightly lower than during 
the subsequent collapsing phase. In these studies entropy refers to that of the content of 
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the universe [T6J [T7] and ignores contributions from (quantum) gravity. 

There are, however, important problems with these models, including the lack of treat- 
ment of singularities and the un-corroborated assumption that the bounces themselves leave 
the entropy of the universe unchanged. The consideration of oscillatory models within a 
quantum cosmological framework, on the other hand, not only allows singularities to be 
avoided, but also introduces many more quantum degrees of freedom, thus allowing the 
question of entropy to be considered in a different light. 

This is the setting we consider in this paper. We will analyze the recollapse in detail, 
which is a semiclassical regime but, crucially, still described in terms of a quantum state. 
We especially focus on the evolution of state parameters through the recollapse, which 
provides insights to the question of what their generic change may be. In particular, we 
are interested in how strongly fluctuations of a generic state respect time-reversal symmetry 
for time reflections around the recollapse point. If fluctuations are symmetric in this sense, 
there is not much change between the pre- and post-recollapse phases. A violation of 
the symmetry, on the other hand, would provide a measure for the generic change of the 
quantum state in the recollapse phase. The analysis is thus complementary to what has 
already been studied for the bounce [HI [19] : Can the quantum state after the recollapse 
be very different from what it was before? Especially in the presence of many cycles, this 
question is important for understanding the viability of oscillatory cosmological models 
over epochs long compared to the life time of individual cycles. 

For technical reasons, we shall take a free massless scalar as the matter source in all 
models considered in detail here. However, we shall also demonstrate the robustness of our 
claims under the inclusion of potentials. The free scalar has the advantage that it can be 
used as a global internal time parameter and thus gives rise to true Hamiltonian, rather 
than constrained, evolution. Any non-constant potential or even a mass term would spoil 
this feature. (Here we refer to the classical situation. We will later encounter and entertain 
the possibility of genuine quantum variables as a measure for time even in situations 
where no obvious classical clock may exist.) Moreover, in the absence of a cosmological 
constant and for flat, isotropic space, this matter content provides an exactly solvable model 
even after quantization (loop or otherwise) [TT]. Thus, there are no dynamical quantum 
corrections whatsoever in this case; the system is harmonic and presents the simplest and 
most controlled model of quantum cosmology. (There may, however, be quantum geometry 
corrections of kinematical type which give rise to a bounce in loop quantum cosmology. 
But they turn out not to spoil the dynamical solvability However, this exact model 

does not allow a recollapse, and we therefore have to add extra ingredients and with them 
non-trivial quantum corrections. Nevertheless, the resulting systems will be manageable 
and provide key contributions for highly controlled cyclic models. While there is no scalar 
potential in the main part of the paper, we verify that in fact our results remain robust in 
presence of general non-zero potentials. Moreover, our analysis provides a starting point 
to analyze equations in the presence of a potential perturbatively. For the bounce, such 
equations are developed in [131 E], which in some cases even allow conclusions valid to all 
orders in the potential and in quantum moments [20J. Since our main question is about 
limitations to the symmetry of fluctuations around cosmological turning points, a highly 
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controlled model is reliable as any limitation there would only grow if the model becomes 
more complicated. (See also [18j[T9] i n this context.) In Sec. 13.31 we will comment in more 
detail on possible effects of a potential. 



2 Recollapsing models 

We shall confine ourselves to isotropic and homogeneous settings. There are two different 
ways to achieve a recollapsing cosmological model: by including a negative cosmological 
constant or by allowing positive spatial curvature. We shall first describe the general 
scheme of our analysis and then specialize to these two cases. 



2.1 Prescription 

In the presence of a cosmological constant and a free massless scalar field the Friedmann 
equation takes the form 

where p^ is the momentum corresponding to the homogeneous scalar field <fi and can be 
written as 



±« 2 \/ T^Vd 2 + k — Aa 2 
V 47rG 



±2\/^(l ~*)V (2) 



in terms of canonical gravitational variables 

o 2-2x 

V =^n r T ^dP = -f a 2x a with {V,P} = 1. (3) 

The introduction of this pair of canonical variables and the parameters fo and x is motivated 
by loop quantum cosmology and deserves further explanation: While /o will not be of 
much consequence in what follows, we keep it for general reference. It has dimensions 
such that P becomes dimensionless; for x = —1/2, for instance, it has the dimension of 
length and for x = it is itself dimensionless. Its significance lies in the fact that it 
determines a fundamental scale for the loop quantization which becomes relevant at the 
bounce. Moreover, the fundamental length can depend on the evolution of the universe, and 
thus a, if the underlying discreteness of quantum gravity is being refined during evolution. 
This possibility is taken into account by the parameter x, which makes the momentum P 
depend differently on a for different choices of x. How precisely these parameters arise has 
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been discussed, e.g., in [2TJ [22]. The dynamical behaviour of loop quantum cosmology is 
sensitive to their values, but in this paper we will mainly analyze recollapses where effects 
of the loop quantization are not expected to play large roles. We will nevertheless see that 
it is of interest to keep all possibilities, especially of x. For all choices of fo and x, the 
variables used here are canonically related to each other. Nevertheless, some quantitative 
aspects can change, and also equations of motion may be easier to solve for some x than 
others. 

Physically, different values of x correspond to different ways in which an inhomogeneous 
discrete quantum state can be refined during its evolution on microscopic levels [21]. For 
x = the variable P corresponds to an underlying state which has a constant number 
of lattice sites as the universe expands, while for x = —1/2 the state has a constant 
geometrical size at each lattice site and thus requires new sites to be generated during 
expansion. A precise value of x could in principle be determined if one could derive a 
reduced Hamiltonian of an isotropic model from a full, inhomogeneous Hamiltonian (such 
as those introduced in |23j). Since this is not yet available, we have to keep the the value 
of x free and look instead for possible phenomenological constraints. 

Classical solutions as functions of <fi are readily determined from the Hamiltonian 
H oc and its canonical equations of motion in terms of 0, which will be presented be- 
low. Such equations of motion determine the relational dependence of, e.g., V{4>) through 
the Hamiltonian equation of motion dV/d<p = {V,H}. Our main interest, however, is in 
possible effects which may result from the behaviour of quantum states. In particular, 
a quantum system has not only expectation values as free variables, which could be as- 
sociated with the classical variables (V,P), but also fluctuations, correlations and higher 
moments. Dynamically, all these variables couple in a general quantum system. These 
coupled equations of motion can be derived from the usual commutator relations uch as 
d(V)/d(j) = — i/t -1 ([V, H]) or, more compactly, from a quantum Hamiltonian Hq := {H). 
(For details we refer to [2^1 125] or, in the context of cosmological models, [T2l [19] .) Here, 
the expectation value is computed in a state with a general set of moments. As is well 
known, for a general classical Hamiltonian H(V,P) we have (H(V,P)) ^ H((V),(P)) 
where the difference amounts to quantum corrections to the classical dynamics. These 
corrections depend, e.g., on quantum fluctuations or, more generally, on moments 

G a > b = ((V - (V)T(P - (P)) fe ) Wey i, (4) 

of the state used for the expectation values. (In the definition of moments, we assume the 
basic operators to be totally symmetric or Weyl ordered as indicated by the subscript.) 
Upon writing Hq = (H) in terms of expectation values and the moments, we obtain 
the complete quantum Hamiltonian. This in turn generates the Hamiltonian equations of 
motion foiQ V := (V), P := (P) as well as all the moments G a,b . (As before, the equations 
of motion are given by df/d<p = {/,Hq} where {V, P} = 1 and for G a ' b the Poisson 
brackets follow from expectation values of commutators divided by ih.) 

^rom now on, we will mostly be referring to the quantum theory unless stated otherwise, and thus 
drop brackets on expectation values. 
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G PP 


= G°> 2 = 


g vp 


= G 1 ' 1 — 


G vv 


= G 2fi = 



This is the basis for the derivation of effective equations which may provide good 
approximations in regimes where the infinite set of all moments can be truncated to finitely 
many variables. In the following we shall only consider the second order moments which, 
for better clarity, we denote as 

(P 2 ) - P 2 (5) 

-{VP + PV) -VP (6) 

(V 2 ) - V 2 . (7) 

Their Poisson brackets can be then derived as in 

{G VV ,G PP } = {{V 2 )-V 2 ,(P 2 )-P 2 } 

1 IP . IV *. - 41/ P - - 

= 2(VP + PV) -AVP = AG VP . (8) 

Similarly, 

{G VV ,G VP } = 2G VV and {G vp , G pp } = 2G PP . (9) 

Such Poisson brackets, when used in dG a,b /d(f) = {G a,b ,HQ}, determine the evolution of 
the quantum variables of a state. This demonstrates how effective equations are able to 
go well beyond simple corrections to classical equations, which will be made ample use of 
in this article. 



2.2 Negative cosmological constant 

For A < 0, k = 0, our system has the classical Hamiltonian 

H = (I - x)V^P 2 + |A|/ 2 (87r 7 G(I - x)f V/3Y 1+2x y^- x \ (10) 

for 0-evolution, i.e. = 27a/87tG /3H (a specific sign has been chosen here for the square 
root; the other choice simply amounts to replacing with —(f)). The factor in can be 
eliminated by redefining 0. Evolution is analyzed best for x — —1/2, in which case 

H = \vjF^Wl (If) 

is linear in V . The corresponding quantum Hamiltonian, including moments of second 
order, is 

which includes the quantum moments G pp , G VP in correction terms. Higher moments are 
ignored here, and G vv does not occur thanks to the linearity of H in V . (For A = we 
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have the solvable free system, in which no coupling terms between expectation values and 
moments arise [H].) The quantum Hamiltonian determines the Hamiltonian equations of 
motion 

^ = 3 VP -W VP C pp \\\\f* (13) 
d<P 2^pT^\K\JI 4 |A|/o (P2 + | A | / 2 )5 / 2 ^ + 2 |A|Jo (P2 + | A | /o2) 3/ 2 W 

% = -fv^^-!i A ^(^SF- (14) 

Quantum fluctuations appear here in coupling terms and are themselves dynamical, 
subject to equations of motion 

Ar pp P 

—j— — = -3 - = G pp (15) 

d0 V^ + IAI/o 2 

dG VP 3 , . , „ 1/ 



5l A l/o , D2 I IAI , 2 vs/ 2 gPF ( 16 ) 
3 l A l/o 77^ ■ hrPK^ GVP + 3 7 ==^= G yv . (17) 



d0 2' '^(P2 + | A | /2) 3/ 2 

# MJo (P 2 + |A|/o 2 ) 3/2 ^ ' V^ 2 + |A|/o 2 

These equations satisfy 

d ( G yy G PP _ {G v P) 2) = o 

such that a state initially saturating the (generalized) uncertainty relation 

g vv g pp _ (g vp )2 > ^ (lg) 

will keep saturating it. Such a state would be considered a dynamical coherent state whose 
properties can be analyzed by our equations. In what follows, however, we will not restrict 
states to be on the saturation surface although they certainly must satisfy the uncertainty 
relation. 

If we first ignore all moments and their quantum back-reaction, we find the classical 
solutions 



Piassicai(0) = Pq cosh(3(0 — <f>o)/2) + y Pq — |A|/q sinh(3(0 — 0o)/2) (19) 

Klassical(0) = V ^ y + BI . (20) 

°-P o sinh(3(0 - o )/2) + y/Pg + |A|^cosh(3(0 - fa)/2) 

The volume has a turning point, and we can simplify expressions without loss of generality 
by choosing our initial values there, i.e. P = P(0o) = and shift such that O — 0. 
Then, we have simply 



Pclassical(0) = ~ V l A l/o Sinh(30/2) (21) 

. . Vn 



These solutions describe the recollapse of a universe with a past and a future singularity. 
Analytical solutions of equations amended by quantum geometry effects, where the singu- 
larities are replaced by bounces and thus provide cyclic solutions, have been derived e.g. 
in [26]. However, quantum back- reaction effects, which complicate the analysis, were not 
included in the equations used there. 

In a next step, we can solve the equations of motion (|T5|) - (fT7j) approximately by assum- 
ing the classical solutions for P and V. Thus, we are still ignoring quantum back-reaction 
effects at this stage, which if present would imply that the moments back-react by the 
coupling terms in ( fl3l) and ( fl4l) and change the classical solutions. For small fluctuations, 
this will be a good approximation, and solutions obtained for the moments will allow us 
to check self-consistently for how long in it will remain valid. 

It is then easy to solve for G pp , to give 

G pp (<f)) = G pp cosh 2 (30/2), (23) 

which shows that G pp is inversely proportional to the volume squared, and which in turn 
allows to solve for 

G^(0) = GV + ^ogT ^(30/2) (24) 
VW/o cosh(30/2) 1 ; 

With this, one can finally solve for 

G vv + 2 tanh ( 3 0/2) + tanh 2 (30/2) 

qw vi A i/° ^ / 25 \ 

cosh 2 (30/2) ' 1 ' 

With quantum back-reaction to second order in moments, i.e. solving the full equations 
(!T3l-(fT7l) without starting with the classical solutions, the equations are more highly cou- 
pled. One can derive some solutions by dividing iHM by ( fl5l) . thus providing a differential 
equation for P{G PP ): 

dP P 2 + |A|/ 2 , 1 |A|/ 2 



+ 7 o,o 7,,^ - (26) 



dG pp 2PG PP 4P(P 2 + |A|/ 2 ) ' 
This can be written in a simpler form thus 

d(P 2 + |A|/ 2 ) „, , IAI , 2 , 1 A ., r 



dlogG pp 1 uu 2P 2 + |A|/ ( 

whose solution yields 



^ 2 + |A|/o 2 +o o 2 I (27) 



P = V-|A|/o 2 + V4G pp y - \A\G pp (21 

such that 



'P 2 + |A|/ 2 = Hc(G PP Y - |A|/ 2 G^, (29) 
with a constant of integration c. 
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2.3 Positive spatial curvature 



With A = but k — 1, the system is simplest to solve for x — 0, which makes it again linear 
in V . The quantum Hamiltonian is then the same as before, ( fl2|) . with A replaced by — 1 
(and a missing factor of 3/2 arising from 1 — x in the Hamiltonian, which simply rescales 
<ft). We can thus immediately take over the solutions already found. For other values of 
x, the equations are more highly coupled and do not allow simple solutions. Nevertheless, 
we can use the solutions already provided to find information also about these systems by 
simply replacing (V, P) in the x = 0-solutions by 

V-:=4«l-*)W-> , ^= ((1 _ x)G ; ^/(-.) = f /( G ^. (30) 

(We have chosen the factors of G and fo such that P has the same dimensions as fo, which 
will be useful later.) This has to be done also in the moments, i.e. we will obtain their 
solutions not for G vv , say, but for G vv . These are not directly the fluctuations of our 
basic variables for x ^ but they still give important information about the spreading and 
other properties of states. For instance, we will determine the correlation G VP instead of 
G VP . Both parameters contain equally interesting information about squeezing and the 
symmetry of fluctuations around the recollapse. In particular, if G vv is not symmetric 
around the recollapse, then nor will be G vv . 



3 Implications 

Several conclusions can be drawn from the solutions found to the given order. 



3.1 Volume ratio between recollapse and high curvature regimes 

Our solutions correspond to state parameters in a Wheeler-DeWitt quantization because 
we use elementary variables (V, P) which are assumed to be quantized to well-defined 
operators. Those operators, together with the Hamiltonian, then determine the dynamics. 
The latter have not been written explicitly here, but they are the central ingredient to 
Hamiltonian equations of motion via the Poisson brackets of quantum variables such as 
© and ®. 

The Wheeler-DeWitt quantization does not easily solve the singularity problem. For 
models without quantum back-reaction effects, i.e. spatially flat models sourced by a free 
massless scalar, (V) simply follows the classical trajectory into the singularity. On the 
other hand, in general models such as those considered here, there are quantum back- 
reaction effects which one may expect to become stronger as the solution for V approaches 
zero — the classical singularity. This could stop V altogether, or delay its approach to 
zero sufficiently strongly such that zero would not be reached in a finite amount of proper 
time (but possibly still finite in <p). However, this is difficult to analyze if all moments are 
required, and unlikely to result in a generic resolution of singularities. 
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A loop quantization does provide a natural solution of the singularity problem in 
isotropic models, but it requires one to use a different set of basic variables. (At a ba- 
sic level, singularities in homogeneous and spherically symmetric models have been shown 
to be absent by allowing general wave functions to be extended through classical singu- 
larities [27J EE! EHl ED]. More specific examples for bouncing wave packets are derived in 
[7JEEI]- For a discussion and comparison of results concerning singularities see [2].) While V 
would still be represented as an operator in the quantization, the curvature (or connection) 
component P is not. Instead, loop quantum gravity is based on a quantum representation 
in which only holonomies of the Ashtekar connection are represented, in this way providing 
the kinematical structures for a well-defined, background independent quantization of full 
gravity [31j [321 [33]. In the cosmological models studied here, this means that it is not P 
which is part of the elementary algebra but exp(i/zP), for arbitrary real fi. (Note it is P 
which enters here, rather than P of ([30]) because x represents the freedom in the refinement 
of a discrete underlying state and thus determines the form of holonomies in a reduced 
isotropic setting [21J. This is, in fact, the main reason why we allow for different values of 
x.) Using the exponential instead of an expression linear in P changes the basic algebra 
as well as the Hamiltonian in particular at large P. In a flat, isotropic model with a free 
scalar field, the classical singularity is then resolved and replaced by a bounce. 

To study the oscillating models we need to consider a combination of bounces and rec- 
ollapses which is more complicated because of the structure of required quantum evolution 
equations. Nevertheless, one can study cyclic solutions by patching together bounce and 
recollapse phases. For small curvatures, we can use the equations and the corresponding 
solutions provided in this paper to an excellent approximation, even for a model of loop 
quantum cosmology. However, we can use this only when P is not too large and have to 
cut off our solutions at the latest when \P\ ~ 1. (At this point, the precise value of fo 
would set the corresponding scale for d.) This leaves only a finite range of sizes for the uni- 
verse between this high curvature regime and the recollapse. The high curvature regimes 
can also be described by effective equations, which are in fact precise without quantum 
back-reaction, but require a different set of basic variables [11] . 

For a negative cosmological constant, we have the ratio Vq/V\p\=i = a/1 + 1/|A|/q . 
Thus, for a small cosmological constant compared to f ~ 2 , the ratio is huge. Since fo 
arises from quantum gravity and has the dimension of length in this case which is based 
on x = —1/2, f should take a value near the Planck length. Thus, |A| must only be 
small compared to a Planckian value which can safely be assumed to be the case. For 
the closed model with x = 0, on the other hand, we have Vq/V\p\ = \ = a/I + 1/ fo with a 
dimensionless fo. In this case, there are no strong reasons to expect quantum gravity to 
provide a value of fo small compared to one (without reference to a second scale larger 
than the Planck length, which should not appear in the basic variables V and P where fo 
enters). This is certainly not enough for a macroscopic universe which has to grow large 
out of the high curvature regime. For this reason we have to use other values for x in this 
case: then, \P\ = 1 is reached at much smaller values for P as provided by our solutions. 
Although the qualitative behaviour is unchanged compared to other x, changes in x have 
an important quantitative implication (which was first emphasized in [10]). For example 
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for x = 



—1/2, the high curvature regime starts at 



cosh(0) ~ - 




108 C + 12 V-12 + 81C 2 + 



y/w&C+ 12^-12 + 81^ 



2 



where C = (GfoVo) 2 ^ 3 //^ . For large Vo, this is approximately cosh(0) « Vq (or, for 

general a; ^ 0, cosh(0) « V^ 1 ^ ). Thus, the ratio Vo/V\p\=i ~ Vq -3 ^ 1-3 ^ is no longer 

constant and grows with Vo for negative x. For x = —1/2, the ratio is given by V 2 ^ 9 which 
is large enough for large Vo, leaving ample room for a growing universe. 

3.2 Quantum back-reaction effects 

From our solutions we can determine whether quantum back-reaction effects are strong 
around the recollapse. As one can easily see, there are no possible divergences in the 
equations of motion f[T3"|) and (TH1) which would enhance the coupling terms. Quantum 
back-reaction effects can only be strong if the quantum variables are large, which can be 
avoided at least for some time by choosing a semiclassical initial state. Thus, the equations 
to the order provided here are reliable to a high degree and can be used to determine the 
state properties around the recollapse. In particular, our equations of motion and solutions 
for quantum variables themselves can be used to see how long the approximation remains 
valid. 

3.3 Evolution of the spread 

Of particular interest is whether fluctuations depend strongly on or remain nearly con- 
stant during the evolution. If they change rapidly, the behaviour of neighbouring cycles 
would be noticeably different from each other because the state would have changed sig- 
nificantly. In scenarios with a large or an infinite number of cycles, large differences should 
even be generic between widely separated cycles. 

As we have seen, G pp is always proportional to the inverse volume squared when quan- 
tum back-reaction effects can be ignored. Thus, curvature fluctuations must be symmetric 
around the recollapse and do not change significantly: At any volume after the recollapse 
we have the same G pp as at the same volume before. For the other quantum variables, 
however, the situation is different. Ignoring products of quantum variables, we can rewrite 
( TT7j) approximately as 




(31) 



for a negative cosmological constant with x 



1/2 or 




(P2 + /o 2)3/2> 



(32) 
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for a closed model with x = and A = 0. This shows that G vv would be a function only 
of V, and thus symmetric around the recollapse, if G VP = 0, i.e. the state is unsqueezed. 
One may assume this as an initial condition, but G VP itself is dynamical and subject to the 
evolution equation ( JT6l) . Its time derivative cannot be zero since, thanks to the uncertainty 
relation, G pp is non-zero unless volume fluctuations diverge. Even an initially unsqueezed 
state will become squeezed after some time, and thus also affect the volume fluctuations. 

Even if G vv /V is not constant, G vv may be symmetric around the recollapse but 
behave differently with respect to V. In fact, (125]) shows that G vv is symmetric around 
the recollapse if Gq P , i.e. the correlation at the recollapse, vanishes even though G VP would 
become non-zero away from the recollapse. But since this happens only under the special 
condition of Gq P = 0, it could generically be satisfied only in one cycle of an oscillatory 
universe. 

From ff24l) . we can estimate the change in squeezing per recollapse by 

lim G VP {cf ) ) - lim G VP (<P) = (33) 
0-+-oa v|A|/o 



as an upper bound. The change may be small for small fluctuations G pp , but is enlarged 
by a factor of V (as well as 1/ ^/\A]f in the presence of A < 0, which is large given that 
|A|/q is small; if the recollapse is triggered by positive spatial curvature, we have the same 
formula with A set to —1). In a large universe, this change can be quite significant. Note 
that in (13"3"|) we have used <fi — > ±oo, and thus a range which includes the high curvature 
regimes where the equations have to be amended by effects of the loop quantization and 
the specific solution would change. We can take this into account by reducing the range 
of 0; however, this does not change the result but only affects the numerical factor in 
the change of squeezing. There is thus a significant change during the classical recollapse, 
irrespective of how the high curvature regime is dealt with. For instance, we have 



t/ r<PP 

|sinh(30/2)=l - <^0 + —fKffu ' ^ 

^| A I Jo 



whose numerical coefficient is different, but which still carries the large factor of Vq. In fact, 
the tanh-behaviour of G VP demonstrates that the greatest change in correlations occurs 
near the recollapse. 

To quantify the production of squeezing during recollapse phases, it may be helpful 
to transform the solution for G VP ((p) to proper time rather than using the relational for- 
mulation with respect to (p. The relation between proper time r and <p can in general be 
complicated, but can easily be obtained for x = —1/2 by integrating 

- = - = -cosh(30/2) 

to obtain 

2 

</>(t) = -arsinh(tan(3p0(r - t\)/2V q )) . 
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Without loss of generality, we chose to vanish at T\ 1 which may be different from the 
recollapse time Tq. The whole range — oo < <p < oo corresponds to a finite proper time 
interval — Vo7r/3p0 < r — t\ < Von/3p^. This highlights the fact that we are not including 
effects of the loop quantization, such that the endpoints of the 0-range, where the volume 
vanishes, correspond to future and past singularities a finite proper time away. 
Inserting this in the solution (T24j) for G VP , we obtain 



which shows the growth of squeezing in proper time during each recollapse (which is in 
fact monotonic in the given range of r). 

Starting with an initially unsqueezed state it may seem that for many cycles the state 
remains almost unchanged from cycle to cycle. Its volume fluctuations may always seem 
to attain nearly the same size at the same volume. However, this is so only because of the 
special initial state chosen, from which squeezing builds up slowly. For small G , (l3Tj) 
and (1521) . respectively, show that G VP /V is nearly constant in both cases considered. The 
change in volume fluctuations before compared to after the recollapse seems insignificant 
from cycle to cycle but becomes noticeable over many cycles. Moreover, if the initial state 
had already had some squeezing, volume fluctuations relative to volume would change much 
more rapidly. In this way, the choice of initial state can strongly influence the long-term 
behaviour. 

In a cyclic model, it is especially important to ask what significance one should attribute 
to the choice of initial state. Is it to be posed in "our" cycle, and if not, how many cycles 
ago? If we could have observational input on properties of the state, we could certainly 
pose an initial condition in our cycle and see how the state evolves to or from there. 
However, state properties are hardly under control, and this possibility remains elusive. 
We thus have to pose initial conditions many cycles ago based on some general principle 
of emergence, but we never know how many. Thus, even though we know that an initially 
unsqueezed state builds up squeezing only slowly, this does not say much about the present 
state if we do not know how many cycles ago the state was unsqueezed. 

An interesting question is whether in a cyclic model one generically expects to have 
a finite or an infinite number of past cycles. The problem with the finite case is that it 
does not resolve the origin question. In the emergent scenarios [3U ESI EU E7J EE], as well 
as some other such models, the universe is assumed to have undergone an infinite number 
of past cycles so as to remove the question of the origin. In that case any given cycle 
would have an infinite number of precursors and generically we therefore have to expect 
the current state to be squeezed. (We will argue in the next subsection that bounces do not 
affect the qualitative behaviour of the squeezing, especially its monotonicity) . The question 
then is how the squeezing in a generic cycle is determined. If each cycle produces the same 
amount of squeezing, a generic cycle would have infinitely squeezed states, which could not 
be semiclassical. However, as (1331) shows, the amount of new squeezing per cycle depends 
on the recollapse volume Vq of that cycle. For growing cycles, as in the emergent scenario, 




(35) 
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the change in squeezing is initially small and approaches zero for cycles in the infinitely 
distant past. Depending on the precise scenario, the sum of all squeezing contributions 
may converge, such that a finite value results for a generic cycle. Whether this is the case 
and what this precise value could be depends on which concrete model one is using, and we 
will not follow this route here. It is, however, interesting that this in principle allows one 
to restrict the possibilities for emergent scenarios by the amount of squeezing they would 
predict. 

Another interesting and related question is that in the emergent models the eventual 
non- uniformity of cycles is produced by a non-constant potential. (In initial regions where 
the potential is flat the universe would just periodically oscillate around the center point; 
the eventual asymmetric emergence is induced by a non-trivial change in the underlying 
potential.) This raises the question of what happens to these models when treated quan- 
tum mechanically. Taking the case of a negative cosmological constant, corresponding to 
a constant negative potential, as a guide suggests that, even though in the flat regions 
of the potential there is a classical symmetry expressed by the exact periodicity in the 
dynamics, we nevertheless acquire a quantum mechanical asymmetry due to the evolution 
in squeezing. A more complicated question is what happens in the regions where there is 
already a classical asymmetry induced by the non-flat potential. 

To be specific, let us look at the closed model with x — 0, while including a scalar 
potential W(<f>). In this case, will no longer serve as a global internal time, but it is still 
a good indicator of local internal time in phases where <f> is monotonic (i.e. outside zeros 
of p^). In this way, we can still draw conclusions for the behaviour of quantum variables 
near a recollapse. In this case we have the Hamiltonian 

H = V^JP* + fi - 87r 7 G^(0)/ o 3 \//3, (36) 
and a corresponding quantum Hamiltonian 

= V V P +/o -^GW(<f>)f V/3 + - 2{P , + n _^ lGW{mv/?>)m G 
, P(P 2 + / 2 - A^GW^V/3) QVP 

G' 



{P 2 + f 2 _ s^GW^flV/^ 
4tt 7 GH/(0)/ 3 (P 2 + / 2 - 2-K 1 GW{<P)f$V)/3 nVV 



(p2 + /2 _ 87T 7 GW(0)/ o 3 l73)3/2 



expanded to moments of second order. In contrast to the previous cases, this includes not 
only the quantum moments G pp , G VP but also G vv in correction terms. The quantum 
Hamiltonian then determines equations of motion, which for G VP results in 



dG VP V(f* - 87r 7 GW(0)/ o 3 y/3) 



:G 



pp 



d0 " (P 2 + / 2 - 87r 7 GH/(0)/ 3 \//3) 3/2 

™*™Vv. (38) 



8tt 7 GH/(0)/ 3 (P 2 + /o 2 - 27r 7 GH/(0)/ 3 ^)/3^ 



(P 2 + fi - 87T 7 GH/(0)/ 3 K/3) 3 / 2 
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Since P 2 + /g — Sn-yGW ((p) f$ V/3 is required to be positive and P is small near a recollapse, 
the sign of this expression remains unchanged compared to the free model. Thus, inclusion 
of a potential does not change our monotonicity result. Notice that we have not assumed 
the potential to be small since the analysis involves only an expansion in moments rather 
than in W ((/>). The rate of change of correlations depends on the value of the potential, 
but it has a definite sign: G VP is either growing or decreasing during a recollapse phase. 
A varying potential will affect the rate by which G VP changes, and thus lead to different 
absolute changes in squeezing before and after the recollapse. But correlations will always 
change, and thus our qualitative discussion remains unchanged in this case. 

In the cyclic models with many cycles one can only draw conclusions from the con- 
sideration of generic rather than special initial states. Thus one needs to consider the 
consequences of generic initially squeezed states, rather than special unsqueezed initial 
states. While a state may have been uncorrelated at some time, we cannot know how 
many cycles ago this may have been, or after how many cycles it may be so in the future. 
For statements relevant to a single cycle, which are the only ones with a chance of being 
observable, it is not legitimate to use special initial states which are known to change be- 
tween cycles. In fact as can be seen from Eqs. fl3T|) or fl32l) . there is no strong bound on 
the change of volume fluctuations relative to volume from one cycle to the next without a 
sharp limit on correlations. Quantum properties of the collapse phase can thus differ from 
those of the expansion phase. As (1251) shows, the time-asymmetric term has a single factor 
of Vb, while the last term is multiplied with Vq . One can thus expect that the asymmetry is 
not pronounced strongly for a universe of large recollapse size V , but the precise behaviour 
depends also on the moments. Then, the last term containing Vq is suppressed by a factor 
of G pp which must be small near the recollapse where P = 0; see Fig. [1] for a numerical 
example. Moreover, over several cycles the change in quantum properties will add up. 

Correlations in a semiclassical state are bounded, and so G VP is restricted but may 
certainly vary. And as long as it can easily be non-zero and affect the behaviour of single 
cycles, it must be taken into account in cyclic models with many cycles. Moreover, in 
addition to the recollapse phases discussed here, squeezing has a similar influence on the 
asymmetry of fluctuations around the bounce [T2]. This shows that the different cycles of 
a universe can indeed be very different from each other, even though they are connected 
by deterministic evolution of an underlying state. The generic behaviour of quantum 
properties is much more subtle than the assumption of unsqueezed states would suggest. 
Current knowledge is insufficient to determine what came before, or what will come after. 

3.4 Entropy 

A central question in cosmology is how to successfully define a notion of cosmological 
entropy, and a number of attempts have been made in this direction. This is in turn hoped 
to provide a notion of cosmological arrow of time. The notion of entropy is connected to 
that of information associated with the degrees of freedom considered. In addition to the 
usual thermodynamic entropy which is normally associated with matter/energy degrees 
of freedom of the constituent components of the universe [17J, possible notions of entropy 
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Figure 1: An example of a recollapsing universe which grows to large volume. Plotted 
are the volume expectation value (V)(<j)) as well as the fluctuations around it. While 
the detailed behaviour of the fluctuations cannot be discerned from this total plot, the 
asymmetry around the recollapse is clearly visible in the zoom shown in Fig. [2j 
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Figure 2: The recollapse phase of Fig. [T] in more detail. The central line is the volume 
expectation value, and solid lines around it illustrate the spread AV of a state which is 
symmetric around the bounce. Dashed lines show how asymmetric the volume fluctuations 
can be if the state is correlated at the recollapse. Initial conditions are set at the recollapse, 
where in units with h = 0.2 we have V = 10 4 , G% v = 200, G pp = 10" 4 and G% p = 
for the solid lines and Gq P = 0.1 for the dashed lines. The latter state thus saturates the 
generalized uncertainty relation (fl8|) . 
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associated with the geometrical [39] and gravitational [ID] degrees of freedom of the universe 
have also been put forward. 

Motivated by the thermodynamical notion of entropy and the associated second law 
of thermodynamics a necessary, but not sufficient, condition that has often been required 
of general notions of entropy is that of monotonicity in time. An important step has 
therefore often been to look for variables defined in terms of the underlying dynamics 
that evolve monotonically. In addition to notions of entropy associated with classical 
degrees of freedom, one would also expect entropic measures associated with the quantum 
mechanical degrees of freedom. An immediate question that any such measure needs to 
answer concerns the nature of its relationship with the thermodynamic measure of entropy. 
In particular an important question in the case of recollapsing/oscillating cosmological 
models is: do (or should) the expanding and recollapsing evolutionary phases possess 
oscillating or monotonic entropies? Furthermore, how should the entropy associated with 
different cycles evolve in oscillating models? 

This question has in fact been the subject of a long standing and intense debate con- 
cerning the relation between the so called thermodynamical and cosmological arrows of 
time [HI [39] . The question is whether the observed asymmetric (monotonic) thermody- 
namical time arrow in the current expanding phase of the universe has a counterpart in 
cosmology, particularly in a recollapsing universe. A number of studies have been made in 
this connection [121 H31 01] . Given the absence of a dynamical explanation for the observed 
asymmetry in the universe, most such studies assume that the observed thermodynamic 
arrow of time must arise from the boundary conditions of the universe |44j. 

Our results above seem to indicate that the degree of the squeeze of the quantum gravity 
state may provide a notion of entropy purely associated with quantum degrees of freedom 
To the best of our knowledge, this is a new possibility not considered before. (Relating 
entropy to the squeezing of a matter state, however, has been considered in the context of 
particle production; see e.g. [121 SSJ, 07J 0H SSI EQJ-) As can be seen from ffl6|) and fl38l . 
the squeezing of a state is strictly monotonic in time during expansion, recollapse and 
contraction of a cycle in the models considered. This demonstrates that even in isotropic 
models, which include the microscopic dynamics only in a highly averaged form, quantum 
aspects prevent one from viewing a collapsing universe simply as a time-reverse of its 
expansion. The quantum theory's arrow of time cannot reverse at the recollapse. 

Unfortunately, it is difficult to follow its evolution through a bounce because this phase 
can only be described in a different set of basic variables ( J = V exp(iP) for P) which make 
the equations solvable. For classical variables, these are easily translated into each other. 
But the transformation is non-linear, such that moments transform in a highly complicated 
way. In any case, the change in squeezing is nevertheless generic because it is unlikely that 
the bounce will restore fluctuations to precisely the value of the preceding cycle. Moreover, 
one can roughly estimate the squeezing as it evolves through the bounce. In the bounce 
phase, only operators such as J := V^exp(zP) exist and give rise to a solvable evolution. 
Moments between V and J can thus be computed exactly [12] , but it is difficult to transform 
between the V-J and V-P moments. However, the bounce happens near P w 7r/2, and 
with 5P := P — tc/2 we have, up to reordering, Re J = (V cos(P)) = (Ucos(7r/2 + SP)) = 
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-(VsinSP) (VSP) = -(V(P - 7r/2)). Thus, ReJ + VP — |V G yp provides an 

estimate for the V-P squeezing as it evolves through the bounce in terms of expectation 
values. Since expectation values are symmetric around the bounce in the absence of a 
potential, not much additional squeezing is generated around the bounces. 

Most of the squeezing is thus generated in the recollapse phases, which resembles recent 
results for cyclic models with bounces based on the Hagedorn phase of string theory [5T] . 
In the present context with a quantum measure for entropy in the form of squeezing, this 
may seem counterintuitive given that the recollapse is a much more classical phase than the 
bounce. However, the production of correlations is not so much a matter of quantum versus 
semiclassical behaviour but rather of the dynamics in a given regime. A state may remain 
semiclassical to an excellent degree, and yet receive a significant amount of squeezing. 
Whether or not this happens depends on the equation of motion for G VP , or the underlying 
Hamiltonian. The analysis presented in this article unambiguously shows the production of 
correlations in a recollapse phase even though it is semiclassical. Although our qualitative 
estimates for the bounce phases are difficult to make precise, the monotonic behaviour of 
correlations at small curvature appears to be an interesting and reliable property. 

The precise amount of squeezing depends on initial conditions. If all moments could 
initially be zero, they would remain so and no squeezing would develop. However, this 
initial condition is impossible because the moments are subject to the uncertainty relation 
( {TBI) . Thus, unless the volume uncertainty diverges, G pp cannot be zero in (fIBT) and an 
initially unsqueezed state inevitably develops squeezing over time which can grow large 
over many cycles. It is thus quantum uncertainty, together with the specific dynamics of 
the system, which prevents the existence of perfectly symmetric states. 

There is a sense in which small squeezing presents a special state with a distinguished 
discrete symmetry. Under time reversal, we map <p — 0, P i— > — P and G VP i— > —G VP 
while the other variables remain unchanged. Thus, a time reversible solution would have 
vanishing squeezing which one may view as a special state analogous to low entropy. As 
fl25|) shows, this is obtained for vanishing correlations at the recollapse. However, since 
G VP would generically be non-zero at a recollapse, especially in a cyclic model, there is 
no solution which is exactly time reversible. Again, it is the uncertainty relation as an 
additional condition, which eliminates those initial values which would correspond to time 
reversal solutions. 

4 Conclusions 

We have studied the evolution of recollapsing models within an isotropic and homogeneous 
quantum cosmological framework in presence of a scalar field. To allow a recollapse we 
consider, in turn, a negative cosmological constant as well as a positive curvature model. 
We derive the resulting quantum evolution equations to second order in moments of a state 
and study their effects on the recollapsing dynamics of the universe, i.e. the expanding, 
turn around and contracting phases. These effective equations allow us to observe that 
state properties generically change during the recollapse, making quantum fluctuations in 
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the expansion and contraction phases different. At large volumes as they are realized at 
a recollapse, the change is not as noticeable as it can be for states travelling through a 
bounce [TBI [19], but it is significant especially in a cyclic model with several recollapse 
phases. As in the case of the bounce, the asymmetry of fluctuations is controlled by 
quantum correlations which have often been ignored in previous studies. 

The specific equations analyzed here thus allow us to identify correlations as a quantum 
measure for the change of fluctuations. More precisely, we find that the squeezing of an 
initial state is strictly monotonic in time throughout these three phases for the models 
considered. Importantly, we have shown this finding to be robust under the inclusion 
of a matter potential. Combining these results with the corresponding ones concerning a 
bounce in loop quantum cosmology we have shown that squeezing of an initial state evolves 
monotonically throughout a whole cycle. The absence of perfectly symmetric states is a 
combined consequence of the specific dynamics of the quantum system together with the 
presence of quantum uncertainty. 

Such monotonicity is of potential importance in two regards. Firstly, it sheds new light 
on a long standing intensive debate concerning the (a) symmetry between the expanding 
and contracting phases in a recollapsing universe. As shown here, the contracting phase 
cannot be a time reverse of the expanding phase. Secondly, it motivates the adoption of 
the degree of squeezing as an alternative measure of (quantum) entropy. 

Qualitatively, we also consider the evolution of the squeezing of an initial state in 
emergent nonsingular oscillating universes in which the universe is assumed to have un- 
dergone a large (possibly infinite) number of past cycles. We argue that the consideration 
of the amount of squeezing in the universe can in principle provide some constraints on 
the viability of such emergent models. In any case, given that a generic cycle does have 
non-vanishing correlations, squeezings of the quantum gravity state must be taken into 
account in order to draw reliable conclusions about cyclic models. 
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